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Abstract 

Employing an extended three flavor version of the NJL model we discuss in detail the phase 
diagram of quark matter. The presence of quark as well as of diquark condensates gives rise 
to a rich structure of the phase diagram. We study in detail the chiral phase transition and 
the color superconductivity as well as color flavor locking as a function of the temperature and 
chemical potentials of the system. 

Introduction 

At low temperatures and densities all quarks are confined into hadrons. In this phase the chiral 
symmetry is spontaneously broken by the quark condensates. Raising the temperature, one expects 
that the chiral symmetry becomes restored and the quarks are free. This state is called quark gluon 
plasma (QGP). In the QGP all symmetries of the QCD Lagrangian are restored. For QCD at low 
temperatures and high densities one expects a phase where the quarks are in a superconducting 
state Jl], ||, ||, 0|. All this different phases define the phase diagram of QCD §] in the plane of 
temperature and density. This phase diagram is not directly accessible. QCD calculations are 
only possible on a lattice at zero baryon density. In order to explore the finite temperature and 
density region, one has to rely on effective models. Two types of such effective models has been 
advanced to study the high density low temperature section The first type of models includes 
weak coupling QCD calculations, including the gluon propagators |^|. The second type include 
instanton Ji], 0, || as well as Nambu and Jona-Lasinio (NJL) models j2|, ^|. These models shows a 
color superconducting phase at high density and low temperature. In this phase the SUc(3) color 
symmetry of QCD breaks down to an SUc(2). Including a third flavor, another phase occurs, the 
color-flavor- locked state (CFL) of quark matter |p| pj)| |pf . 

The two flavor results of the instanton approach are reproduced by the model of Nambu and 
Jona-Lasinio (NJL) [jl3||(l5] if one includes an appropriate interaction as was shown by Schwarz 
et al. fl5|| . This model has been extended by Langfeld and Rho |l6| who included all possible 
interaction channels and discovered an even richer phase structure of the QCD phase diagram, 
including a phase where Lorentz symmetry is spontaneously broken. 

The choice of the NJL model is motivated by the fact that this model displays the same symme- 
tries as QCD and that it describes correctly the spontaneous breakdown of chiral symmetry in the 



vacuum and its restoration at high temperature and density. In addition, the NJL model has been 
successfully used to describe the meson spectra and thus is able to reproduce the low temperature 
and low density phenomena of QCD [jlTj jl9| . Thus, this is a model which starts out from the 
opposite direction as the instanton model which is motivated as high density approximation of 
QCD. Therefore it is interesting to see whether the NJL model is able to describe the other phases, 
the color superconducting phase and the color flavor locking observed in the instanton approach. 
The shortcoming of the NJL model is the fact that it does not describe confinement, or more 
generally any gauge dynamics at all. Here we will evaluate the thermodynamical properties of the 
quarks in the NJL model at finite temperature and density and we will discuss the symmetries of 
the different phases. We present numerical results for the calculation of the different condensates. 
For our study of the phase diagram we use one specific set of parameters. We treat the three flavor 
version of the model, including an interaction in the quark-antiquark channel, a t'Hooft interaction 
and an interaction in the diquark channel. We restrict ourself to the scalar/ pseudoscalar sector of 
these interactions. 

The paper is organized as follows: In chapter [j] we will briefly review the NJL model and present 
the Lagrangian we will use. In chapter || we study the quark condensate and the restoration of 
chiral symmetry. In chapter || we add the interaction in the diquark channel and present the 
numerical results for the color superconducting sector. We will have a complete evaluation of the 
phase diagram of the NJL model, including chiral and superconducting phase transition at finite 
temperature and (strange and light quark) density. In chapter^ we present our conclusions. 

1 The model 

The model we use is an extended version of the NJL model, including an interaction in the diquark 
channel. In fact, the NJL model can be shown to be the simplest low energy approximation of QCD. 
It describes the interaction between two quark currents as a point-like exchange of a perturbative 
gluon [^0|||2l]]. Applying an appropriate Fierz-transformation to this interaction, the Lagrangian 
separates into two pieces: a color singlet interaction between a quark and an antiquark (Jzf^)) and 
a color antitriplet interaction between two quarks J£( qq )- The color singlet channel is attractive 
in the scalar and pseudoscalar sector and repulsive in the vector and pseudovector channel. The 
Lagrangian in the diquark sector has two parts, both attractive: a flavor antisymmetric and flavor 
symmetric channel. The former includes Lorentz scalar, pseudoscalar and vector interactions, the 
latter a pseudo scalar interaction only. 

The coupling constants of these different channels are related to each other by the Fierz trans- 
formation. Due to the extreme simplification of the gluon propagator in this approximation, the 
resulting model cannot reproduce confinement which is described by the infrared behavior of the 
gluon propagator. 

The resulting Lagrangian has a global axial symmetry J7a(1), and an extra term ^£a in the 
form of the t'Hooft determinant is added in order to break explicitly this symmetry. The resulting 
Lagrangian has then the general form: 

J? = -S?0 + -^(qq) + -^{qq) + -^A (1) 

where ,£?o is the free kinetic part. 

The interaction part of the Lagrangian has a global color, flavor and chiral symmetry. Chiral 



symmetry is explicitly broken by non zero current quark masses, flavor symmetry by a mass 
difference between the flavors. 

The different interaction channels of this Lagrangian give rise to a very rich structure of the 
phase diagram, which was completely evaluated in the two flavor case by Langfeld and Rho | jl6| . 
Here we will concentrate on the three flavor case. The evaluation of the complete phase structure in 
the three flavor case is a quite difficult task and we will concentrate here on the Lorentz scalar and 
pseudoscalar interactions. In the mesonic channel this interaction is responsible for the appearance 
of a quark condensate and for the spontaneous breakdown of the chiral symmetry. In the diquark 
channel it gives rise to a diquark condensate which can be identified with a superconducting gap. 

Describing the quark fields by the Dirac-spinors q, the Lagrangian we will use here has the 
form: 
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The first term is the free kinetic part, including the flavor dependent current quark masses too/ 
which break explicitly the chiral symmetry of the Lagrangian. The second part is the scalar/ 
pseudoscalar interaction in the mesonic channel, it is diagonal in color. The matrices Xf act in the 
flavor space. The third part describes the interaction in the scalar/ pseudoscalar diquark channel. 
The charge conjugated quark fields are denoted by q c = Cq T and the color (a,/3,7) and flavor 
(i,j,k) indices are displayed explicitly. We note that due to the charge conjugation operation the 
product qij5q c is a Lorentz scalar. This interaction is antisymmetric in flavor and color, expressed 
by the completely antisymmetric tensor e w . Finally we add the six point interaction in the form 
of the t'Hooft determinant which breaks explicitly the Ua{1) symmetry of the Lagrangian. The 
det runs over the flavor degrees of freedom, consequently the flavors become connected. 

The NJL model is non renormalizable, thus it is not defined until a regularization procedure has 
been specified. As we are interested in the thermodynamical properties of the model, calculated 
with help of the thermodynamical potential, we will use a three dimensional cut-off in momentum 
space. This cut-off limits the validity of the model to momenta well below the cut-off. 

The model contains six parameters: The current mass of the light and strange quarks, the 
coupling constants Go and Gs and the momentum cut-off A which are fixed by physical observables: 
the pion and kaon mass, the pion decay constant, the mass difference between r\ and rj , once the 
mass of the light quarks was fixed, as well as by the vacuum value of the condensate < qq > 1 / 3 = 
— 230MeV. The last parameter is the coupling constant in the diquark channel Gdiq- For the 
mesonic sector we will use the parameters of (2^]: a current light quark mass mo, = 5.5 MeV, 
a current strange quark mass too s = 140.7 MeV, a three dimensional ultraviolet cut-off A = 
620 MeV, a scalar coupling constant Gs = 1.835/A 2 and a determinant coupling Gd = 12.36/A 5 .. 
This parameter set results in effective vacuum quark masses of m q — 367.6 MeV , m s — 549.5 MeV 



and the quark condensates are ((qq)) = (-242MeV)~ 3 and ((ss)) = (-258 MeV)~ 3 . 
We perform our calculations in the mean field approach for an operator product 

hfa « Pi((p 2 » + ((pi»/5 2 - «pi))«/J2» (3) 

where ((/?)) is the thermodynamical average of the operator and the fluctuations around this mean 
value are supposed to be small. We will apply this approximation to the products of quark fields 
appearing in the interaction part of the Lagrangian. 

2 Chiral phase transition 

We start our study with an investigation of the quark-antiquark sector and the chiral phase tran- 
sition. The diquark sector is subject of the next section. 

The NJL model displays the right features of the chiral symmetry breaking. On the one hand 
we have an explicitly broken chiral symmetry by the inclusion of a small current quark mass. On 
the other hand, the model describes correctly the spontaneous breakdown of chiral symmetry: the 
existence of a quark condensate, responsible for a high effective quark mass and the existence 
of massless (or very light, if the chiral symmetry is explicitly broken) Nambu-Goldstone bosons. 
Lattice QCD calculations show that at a temperature of w 170MeV the chiral symmetry is restored 
(the quark condensates melt at increasing temperature), a result which is reproduced by the NJL 
model P^]|pq|pOtl- As the region of finite density is not accessible to lattice QCD calculations, the 
chiral phase transition at high density is a subject of speculation. The point and the order of the 
chiral phase transition in the temperature - density plane defines the phase diagram. Here we will 
present such a phase diagram for the three flavor NJL model and a specified set of parameters. This 
phase diagram can be viewed as an approximation of the QCD phase diagram, but we have to take 
into account that the NJL model does not describe confinement (we always have a gas of quarks 
and not a gas of hadrons) and that the degrees of freedom are not the same as in QCD (the model 
contains no gluons) . Here we will focus on the thermodynamical properties of the quarks described 
in the color singlet channel of the Lagrangian (||) , this means the thermodynamical properties of 
the quark condensates and masses. 

For the study of the thermodynamical properties of the quark-antiquark sector we will evaluate 
the thermodynamical potential in the mean field approximation. We start out from the Lagrangian 
in the mean field approximation 

Jf MF = q{0 - M)q - 2G s (a 2 + [3 2 + 7 2 ) + 4G D a/3 7 . (4) 

where Mf is the effective quark mass (defined via the quark condensates ((qq))) 

M f = m of -AGs((q f qf))+2G D ((q fl q fl ))((q h q h )) i±h±h (5) 

= niof + 8m f 

and the quark condensates are written in a short hand notation 



a=((uu)) (3 = ((dd)) 7 ={{ss)). 



(6) 



The mean field Hamiltonian 



H MF = Jd 3 x Yl Wl° d m + 2G s(" 2 + I 32 + 7 2 ) - 4G D a/3 7 ] (7) 



f={uM,s} 

is transformed into an operator H in second quantization using 

d 3 p 



At the moment, the quark condensates are unknown quantities. In order to evaluate them, we 
calculate the grand-canonical potential 



(9) 

with [i being the chemical potential, (3 the inverse temperature and N the particle number operator: 



N = h(p, s, f, c) - h(p, s, /, c) (10) 

where h(p, s, /, c) ,n(p, s, /, c) are the number operators for particles and antiparticles with momen- 
tum p, spin s, flavor / and color c. These operators are defined via the creation and annihilation 
operators for particles n(p, s, /, c) = at g j. c a^ s j^ c and antiparticles h(p, s, /, c) = bt g ^ J>p, a ,f,c- We 
consider the condensates as parameters with respect to which the potential has to be minimized. 
The appearance of the quark condensates breaks spontaneously the chiral symmetry of the original 
Lagrangian. 

In second quantization the exponent of the chemical potential reads as follows : 
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-(E P j + n f )h(p,s,f, c)] 

+ [2G s (a 2 +/3 2 + 7 2 )-4G D a/3 7 ] (11) 



where V denotes the volume we have integrated out. The energy E^j = \J M 2 + p 1 depends on 
the flavor of the quarks and their momentum but is independent of color or spin. The evaluation 
of the grand canonical potential in the mean field approximation gives the result: 

= 2G s (a 2 +f3 2 + 1 2 )-4G D af3 1 

E [ A p 2 d P {E PJ + hn[l + e-^f-^] (12) 

+^ln[l + e -0( B */+f/)]j. 

It has to be minimized with respect to the quark condensates : 

dn MF 



0. (13) 
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Figure 1: The mass of strange and light quarks as a function of the temperature (Ihs) and as a 
function of the light quark chemical potential (rhs) 

We obtain three equations, one for each quark condensate 

n c r A p 2 

((q f q f )) = —Mf y / dp^- [1 - f{E pj + M/ ) - f{E pj ~ (14) 

~ JO toft 

where we defined the Fermi function f(x) = (1 + exp(— fix)) -1 . The equations for the quark 
condensates are coupled (see eq.||). For three flavors we have thus three coupled gap equations 
which have to be solved self-consistently. Their solution, displayed in appendix B, enables us to 
calculate the quark condensates and quark masses at finite temperature and chemical potential 
(density). 

We have to take care about the limits of the theory: The regularization cut-off of the theory 
implies that the chemical potential has always to be smaller than this cut-off and that the temper- 
ature has not to be too elevated: The Fermi function will be smoothly extended to high momenta 
and we have to take into account that all states above the cut-off are ignored by the model. 

The condensate is responsible for the spontaneous breakdown of chiral symmetry at low densities 
and temperatures. At high temperature and density the quark condensate drops (it becomes very 
small, or zero - in the case of zero current quark masses) and consequently chiral symmetry is 
restored (up to the current quark masses). Hence the quark condensate is the order parameter 
of the chiral phase transition. The phase transitions we are dealing with are - depending on the 
parameters and of the density respective temperature - of first or second order or of the so called 
cross-over type and we can classify the phase transition by means of this order parameter. The first 
order phase transition is specified by a discontinuity in the order parameter, for the second order 
phase transition the order parameter is continuous but not analytical at the point of the phase 
transition. The third type, the cross-over, is not a phase transition in the proper sense. Here the 
order parameter does not display a non-analytical point but it shows a smooth behavior. 

In a first step we will consider the chiral phase transition as a function of temperature and 
chemical potential of the light quarks, the strange quark density is supposed to be zero. In figure [j], 
lhs, we plot the mass of the light and strange quarks as a function of temperature at zero baryon 
density for the parameters presented above. 

At zero density we observe a smooth cross-over of the chiral phase transition as a function 
of temperature: at low temperature the chiral symmetry is spontaneously broken, with rising 
temperature the quark condensate melts away and the quark masses approach the current mass, at 



least for the light quarks. For the strange quarks we observe a much smoother transition and at the 
highest temperature we can treat in the framework of the NJL model (approximately 230 MeV) 
their mass is still higher than their current mass . This smooth cross-over we observe only for the 
special case of three non-zero current quark masses. 

At zero temperature we observe for our parameter set a first order phase transition. As a 
function of the chemical potential the light quark mass drops suddenly to a value close to the 
current quark mass. The strange quarks change slightly their mass due to the coupling between 
the flavors. For higher values of the chemical potential the strange quark mass is stable. The light 
quark condensate is too small for a change of the strange quark mass. Only a rise of the chemical 
potential of the strange quarks can drop the strange quark mass further, as will be discussed in 
the last part of this section where we present the extension to strange quark matter. 

A first order phase transition is characterized by the existence of metastable phases, the equiv- 
alent of for example oversaturated vapor. These metastable phases are a solution of the gap 
equation, but their thermodynamical potential is larger than for the stable phase. We show this 
in detail in figure g. On the top we display the quark mass (light and strange) , on the bottom the 
density of light quarks and the thermodynamical potential. The stable phases which minimize the 
thermodynamical potential are shown as dark lines, the metastable phases as light lines. 

For the mass of the light quarks we observe the transition from the stable phase at high chemical 
potential to a state whose mass is larger than its chemical potential, this means to zero density. 
Increasing the chemical potential yields a first order phase transition, i.e. the mass of the quarks 
drops suddenly. This abrupt change in the quark masses gives rise to a jump in the density - for 
a constant chemical potential suddenly much more states become accessible. This implies at the 
same time that certain densities does not exist. In our case the normal nuclear matter density is 
just in this region and there are nice explications for this fact [§7]]||. For the interpretation one 
has to remember that we are talking about a quark gas without confinement. Here, nuclear matter 
at normal density one has to consider as a phase which contains dense droplets of quarks in which 
chiral symmetry is restored, surrounded by the vacuum or a very diluted quark gas (which should 
be confined in QCD). The size of these droplets is not given by the theory, but it is not farfetched 
to identify these objects with the nucleons. 

We observe thus for our set of parameters a first order phase transition as a function of the 
chemical potential at zero temperature and cross-over as a function of temperature at zero density. 
Extrapolating now to the plane of finite temperature and chemical potential there must be a point 
where both kinds of phase transition join, the so called tricritical point. In figure || we show 
this phase diagram at finite temperatures and chemical potentials (lhs, at the rhs as a function 
of density). We display as dark lines the transition by the stable state (or the transition line 
for the cross-over) and as light lines the metastable phases. The tricritical point is located at 
a temperature T — 66 MeV and a chemical potential of fi q = 321 MeV which corresponds to a 
density of p q = 1.88po- The location of the tricritical point depends strongly on the choice of the 
cut-off and of the coupling constant J2j| . 

In figure |^ we plot the quark masses (light, lhs and strange, rhs) as a function of the chemical 
potential of light and strange quarks at zero temperature. We can see the influence of the cou- 
pling between the flavors as already discussed for the light quark chemical potential. The strange 
quark mass drops suddenly at high chemical potentials of the strange quarks p, s and low chemical 
potentials for the light quarks fi q . Once the chiral phase transition for the light quarks has taken 
place (at high values of fi q ), the strange quark mass shows a cross-over transition for high fi s . 
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Figure 2: Detailed representation of the first order phase transition as a function of the light quark 
potential at zero temperature. On the top: the light (Ihs) and strange (rhs) quark masses bottom: 
the density (Ihs) and the thermodynamical potential (rhs). The light lines represent the metastable 
region, the dark lines the stable region (minimization of the thermodynamical potential) . 



> 
a> 



a. 
E 



180 ; 

160 : 
140 : 
120 : 
100: 
80 : 
60 : 
40: 
20: 
0^ 



\ <qq>= 

<qq> ^0 \ 



50 100 150 200 250 300 350 400 

chemical potential u q [MeV] 




0.5 1 1.5 2 2.5 3 

density p /p 



Figure 3: Phase diagram for the mass of the light quarks (chiral phase transition) as a function of 
temperature and the light quark chemical potential (Ihs) and density (rhs). The dark lines represent 
the transition by the stable phases, the light lines the transition by metastable phases in case of a 
first order phase transition. 



Figure 4: Light (Ihs) and strange (rhs) quark masses as a function of the chemical potential of light 
(fjL q ) and strange (fi s ) quarks at zero temperature. 



For high values of fi q and /x s , both quark masses have a value close to their current quark mass. 
With increasing temperature, the phase transitions will take place at lower values of the chemical 
potentials. This is very pronounced for the light quarks (see figure |l|) and less for the strange 
quarks which change their mass quite slowly with temperature due to the high current quark mass 
(compare figure 0) . 



3 Color superconductivity 

In this section we will study the diquark channel. We will see that quarks which have opposite 
spin and momenta condense in the scalar channel into diquarks. This resembles superconductivity 
Here we have in addition a complex structure in color and flavor space. In classical 



superconductivity the condensation occurs close to the Fermi surface. In our case we have to take 
into account that quarks with different flavors may have different Fermi-surfaces. Because the 
coupling between the quarks is quite small, the condensation will only occur if the Fermi momenta 
of the two quarks are quite close to each other. 

In order to calculate the properties of the N.IL model in the superconducting sector, we will ap- 
ply the generalized thermodynamical approach of the Hartree-Bogolyubov theory to quark matter 
(see for example J2(|) described by the Lagrangian (||). 

The Lagrangian (§) in the mean field approximation including the diquark sector reads as 
follows: 
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MF 



A fe7 



2G s (a 2 +(i 2 + 7 2 ) + 4G D a/3j 



Qia ' 



A fc 7t 



| A fc 7 |2 



4G 



DIQ 



(15) 



greek indices denote the colors, latin indices the flavors. 
The diquark condensate is defined by: 



A fc T = 2G DIQ i l5 e a ^e^ k {{q v a a l5 e l ' j ' k e a qf, „,)) = i l5 e a ^e^ k A k ~< (16) 

This diquark condensate occurs for all three colors simultaneously. We note that as in classical 
superconductivity the baryon (or particle) number is not conserved. Hence the electromagnetic 
U em (l) symmetry is spontaneously broken and Goldstone bosons appear in the form of Cooper 
pairs. The diquark condensate carries a color and a flavor index. For a given flavor and color the 
condensate is completely antisymmetric in the two other flavors and colors. The condensate A sr 
is created for example by green and blue up and down quarks. 

The diquark condensate is completely antisymmetric in the color degrees of freedom, a property 
which is only shared by three of the eight Gell Mann matrices which generate the SUc(3)- Hence a 
finite diquark condensate breaks down the SUc(3) color symmetry to a SUc(2) if the mass of the 
strange quark is heavy. The same is true for the flavor sector if the three flavors are degenerated in 
mass. For two flavors only the Lagrangian is invariant with respect to a chiral transformation. If 
the diquark condensates coexist for all three flavors, the chiral symmetry is spontaneously broken. 

Due to the product of two antisymmetric tensors the symmetry is even more reduced if all 
three quark flavors form a diquark condensate. In order to see this, we assume first that all three 
colors (for one flavor) are equivalent. Than we can assume without loss of generality that k = 7 in 
equation ( |l6|) and write the tensor product as: 

e iji e *0i = (S ha S ]:P - 5i,p5 j!a ) (17) 

We see that in this case the rotations in color and flavor space are no longer independent but 
locked. Hence the quarks are in a color-flavor locked phase if all three quark flavors participate 
at the formation of the diquark condensates. 5^ a is the unit matrix of SU(3)cxF in which the 
matrices contain as columns the three flavors and as rows the three colors. The Lagrangian is 
therefore invariant under a SU(3)cxF transformation and consequently the SU(3) color and flavor 
symmetries are reduced to an SU(3)cxF symmetry. For more consequences of the appearance of 
the condensate for the symmetries we refer to the literature |Q. Here we will focus on a numerical 
evaluation of the size of the condensates and the phase transitions at finite temperature and density. 

3.1 Thermodynamics 

As before in the case of the chiral phase transition we will evaluate all condensates and the phase 
diagram by the evaluation of the thermodynamical potential. 

We start by writing the Lagrangian in a more symmetric form, following Nambu who developed 
this formalism for the classical superconductivity |^5|. For this purpose we rewrite the Lagrangian 
as a sum of the original Lagrangian and its charge conjugate: 



^ Nambu 

Jz* + ^ G (18) 
Then the Lagrangian can be presented as a matrix: 
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(19) 



where we suppressed the indices for convenience and defined the term 



^cond = -2G s (a 2 +/3 2 + 7 2 ) + 4G D a/3 7 + 



| A fc 7 |5 



In order to calculate the thermodynamical potential in this notation 

SI = -l3Tr[lnexp(-[3(H Nambu - - ^N c ))] 
we need the particle number operator and its charge conjugate: 
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p, s ,f, c 



X I s ?.' 



at — bff ,St 



P, s ,f, c 



(20) 



(21) 



where we suppressed the explicit dependence of the operators on flavor and color degrees of freedom. 

When calculating the Hamiltonian in the mean field approximation, one can see that it is 
possible to separate H — /j,N into two parts, one for the quarks (operators a and fit) and another 
for the antiquarks (operators b et ¥): 



H-liN=(H- fiN)a + (H- M iV) g 



(22) 



These two parts yield the explicit expressions 
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(24) 



We denoted by H con d the expression — V * ^f C ond and used here the discrete summation over the 
momenta. The expressions have a defined structure in flavor and color, the diagonal terms are 
diagonal in flavor and color, the off-diagonal terms (A) are antisymmetric in color and flavor and 
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(25) 



This normalization factor is due to the fact that we deal with products of spinors for different 
species in the off diagonal terms, of course N(p) — 1 when f\ = fa. The explicit form of this 
matrix including all flavor and color indices is displayed in appendix C. 

In order to calculate the thermodynamical potential, we have to diagonalize these expressions. 
This has to be done by means of a Bogolyubov transformation which determines the energies 
of the quasi-particles and the corresponding quasi-particle operators. From the discussion of the 
symmetry of the diquark condensate we expect two quarks of different flavor and color to form a 
diquark condensate whereas one quark of the third flavor is not involved in forming this condensate. 
This has to be seen in the quasi-particle energy and is confirmed if we evaluate explicitly the quasi- 



particle energies as the eigenvalues of the matrices. The diagonalized operators corresponding to 
(H — fiN)a can be expressed in the form: 

3 

(H - fiN)a A = i E ioia\ ia a.Ai a + E' la a Ala a\ la ) 

where i runs over the flavors and a over the colors. a& ia and a^ Aia are annihilation and creation 
operators for the quasi particles. : 
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where 

E~ = E-fi (26) 

E~ = E s -^ s (27) 

Z = ^A\ s N 2 {p) + (E-+Es) 2 (28) 

Y = l(A 2 qq + 4A 2 qs N 2 (p)+E- 2 + Ef) (29) 

X 2 = 1 -(A 2 qq + (8A 2 q8 N 2 (p) + (E-+E-) 2 )(E--E-) 2 

+ 2A 2 qq (4A 2 qs N 2 (p) + (E-+E-)(E--E-))) (30) 

and gi is the degeneracy. For the calculation of the thermodynamical potential it is not necessary to 
know the exact form of the Bogolyubov transformation which relates the quasi particle operators 
Sa with the original quark operators a. The quasi-particles are still fermions, and that is all 
information we need in order to evaluate the sum over the occupied states. It is just necessary to 
assign the right energies to the operators. We evaluate the thermodynamical potential for the case 
of two degenerated light quarks: 
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(31) 



(32) 



This thermodynamical potential contains the (quark and diquark) condensates as parameters. 



In order to evaluate them, we have to minimize 



— = 0. (33) 



<%/3/» 9A flh 

This minimization yields the gap equations for the quark condensates. For the SU(3) case 
the derivation is given in the appendix D. These equations are coupled, we have to solve them 
selfconsistently. The resulting < qq > condensates may be found in appendix E. 



3.2 Results at finite temperature and density 

For this part we decide to take parameter in ref. [Q '-m 0q = 5.96MeV, A = 592.7AfeV, Gs — 
6.92GeV~ 2 , Gdiq/Gs = 3./4. and too, s = 130.7AfeV . We use the relation between the coupling 
constants, ( Gdiq — 3Gs/4 ), given by the Fierz-transformation (see appendix A), close to ]3C| ] 
(0.73). 

The condensates resp. masses at zero temperature as a function of the chemical potential 
fi q = are displayed in figure |[ On the lhs of this figure we show the light and strange quark 
mass, on the rhs the diquarks condensates. First we have to note that quark and diquark conden- 
sates compete with each other as they are formed by the same quarks. Temperature and density 
determine which condensate dominates. 

When the chiral phase transition occurs (the quark condensate disappears), we observe for 
the light quarks that the superconducting phase transition takes place and we have a diquark 
condensate. As the two transitions are related, they are of the same order. The same scenario 
repeats itself for the strange diquark condensate at a higher chemical potential. In the figure we 
display only the solution which is the global minimum of the thermodynamical potential. 

At a quite low chemical potential (the light quarks have a very small mass, the strange quarks 
are heavy) we have only the light diquark condensate, the diquark condensate including strange 
quarks is almost zero as the strange quarks show a strong quark condensate. Only when the 
strange quark condensate drops and the mass of the strange quarks approaches its current mass, the 
strange diquark condensate appears. We have here the coexistence of the light and strange diquark 
condensate, this is the regime where the chiral symmetry is broken again and color and flavor are 
locked. This happens at a quite high chemical potential, the decreasing diquark condensate for 
even higher chemical potentials indicates that we reach the limit of the model: we are too close 
to the cut-off. The phase transitions concerning the strange quarks are quite close to the limits 
of the models if we suppose the current mass of the strange quark of around 140 MeV. We note 
that due to the relatively small difference between the quark masses, both diquark condensates 
have approximately the same value, for the maximum we get A qq rj A qs w (120 AleV) 3 . At 
zero temperature the chiral phase transition (where the quark condensates disappear) and the 
superconducting phase transition (where the diquark condensates appear) are strongly related in 
of our model. This changes at higher temperatures. There the diquark condensates A qq extends 
to smaller values of the chemical potential whereas we need higher densities in order to form 
a A qs diquark condensate. In addition the diquark condensate become smaller with increasing 
temperature. This is shown in figure || where we plot the diquark condensates as a function of 
temperature and chemical potential. For a given chemical potential we observe - as in the classical 
superconductivity - a second order phase transition as a function of the temperature. 

In a next step we consider the diquark condensates in the /i 9 — ^jl s plane. As already men- 
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Figure 5: The light (dark line) and strange (light line) quark masses and the diquark condensates 
A qq (dark line) and A qs (light line) as a function of the chemical potential ji q = fx s at zero 
temperature. 
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Figure 6: Strength of the diquark condensates as a function of the chemical potential fi q = fx a and 
the temperature. We show as color-levels the strength of the condensates A qq (Ihs), A qs (middle) 
and superimpose both as a contour plot (rhs). 
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Figure 7: As a function of the chemical potential fj, q and /i s we show the quark masses on the upper 
row (light quarks Ihs, strange quarks rhs) and the diquark condensates in the lower row (A qq on 
the Ihs and A qs on the rhs) 



tioned, we expect the formation of a diquark condensate only if there are quarks with similar 
Fermi momentum, independent of their mass. Because Gd is zero the disappearance of the quark 
condensates < ss > and < qq > does not depend on the chemical potential of the other species. 
There is one exception the creation of the strange diquark condensates lowers the strange quark 
condensate and increases the light quark mass. In figure [?| we plot the strength of the diquark 
condensates. Because both light quarks have the same chemical potential a diquark condensate 
A qq between the two different flavors occurs whenever the light quark mass is small. 

The strange diquark condensate exists only in a band where the chemical potential of the light 
and strange quarks are approximately equal. The slight deformation of this band is due to the 
different current quark masses. The width of the band is determined by the coupling strength: if 
the coupling in the diquark sector is strong, the quarks can bind and form a condensate even if 
their chemical potentials are quite different. For a small coupling strength, the chemical potentials 
of the two quarks have to be (approximately, in case of different quark masses) equal in order to 
form a diquark condensate. 

Now we should study the feedback of the formation of the diquark condensates on the quark 
condensates (or the mass) . In ^ we display the masses of the light and strange quarks as a function 
of n q and [i s , A qq appears at the chiral phase transition, when the light quark condensate disappears 
and it is formed by the free light quarks. This behavior is almost independent of fi s . Only if A sq 
becomes finite the lack of quarks for the quark condensate increases the light quark mass. The 
behavior of A qs is generic: When the diquark condensate A sq is finite it takes quarks from the 
strange quark condensate lowering the mass of the strange quark. 



4 Conclusions 

In conclusion, we presented the phase diagram of the SU (3) flavor NJL model extended to the 
diquark sector for a set of parameters which reproduces meson masses and coupling constants. 
We find a rich structure of condensates and regions where no condensate exists. The temperature 




Figure 8: As a function of the chemical potential n q and fi s we show the quark masses (light quarks 
Ihs, strange quarks rhs) to pointed out the effect of the diquark condensate on the quark condensate 
at T=l MeV. 

and density dependence of quark and diquark condensates is calculated in mean field approach by 
minimizing the thermodynamical potential. 

The order of chiral phase transition depends on the value of T and \i where the phase transition 
occurs. At zero temperature the phase transition is first order and at zero chemical potential we 
observe a cross over ( due to finite current quark masses). Therefore there exists a tricritical point. 
Normal nuclear matter density does exist only as a mixed phase of a dense quark phase (where 
chiral symmetry is partially restored) and a very diluted quark gas or the vacuum (where chiral 
symmetry is spontaneously broken). Finally we extended the chiral phase transition to the plane 
of finite strange quark density, relevant for the discussion of the diquark condensates. 

Following the idea that the NJL model can be considered as an approximation of the QCD 
Lagrangian we extend the NJL model by including an interaction in the diquark channel. We find 
that this interaction gives rise to a diquark condensation which is responsible for the formation of a 
superconducting gap. This condensation occurs at low temperature and high density. As this gap 
is formed by two quarks of different flavor, their chemical potential has to be close to each other in 
order to allow for this condensation. In flavor SUf{3) this condensate breaks the SUf(3) x SUc(3) 
down to a SUcxf(3) a phenomenon which has already been observed in phase diagrams based on 
instanton Lagrangians and has been dubbed color-flavor locking. We can conclude that two quite 
differently motivated phenomenological approaches to the QCD Lagrangian provide a very similar 
phase structure. 

The diquark condensates do not exist at temperatures expected to be obtained in relativistic 
heavy ion collisions. In neutron stars, which have a high density and a very low temperature they 
could be of relevance. 
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A Fierz Transformation 

Following Ebert (pl|) we have the following relations in color and flavor space: 
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This leads to the following relation between the differents coupling constant: 



qq channel G S ca = ^g 



qq channel G DIQ = -g 



B < qq > condensates 
B.l light quark condensate 
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B.2 strange quark condensate 
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C Matrix 



The total matrix can be separated into 4 submatrices 
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These submatrices are given by: 
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D Thermodynamical potential 

D.l formal derivation of the thermodynamical potential 
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D.2 N(p) derivatives 
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E <qq> condensates 
E.l light diquarks condensate 
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E.2 strange diquarks condensate 
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